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Abstract TUTI 


Inverse equivalent source solutions are powerful approaches for antenna field 
transformations and they provide in particular very detailed diagnostic insight into the 
radiation behaviour of the considered antenna. Commonly such techniques work with 
electric and/or magnetic surface current densities, where the field observation equations are 
sometimes augmented by a null-field condition for the region external to the Huygens 
surface. The surface current densities obtained with such a null-field condition are known 
as Love currents and they are directly linked to the tangential electric and magnetic fields 
on the Huygens surface. Also, it 1s clear that these currents radiate in general (for convex 
surfaces) only away from the antenna under test (AUT) and are, therefore, very suitable for 
antenna diagnostics. As an approximate and inexpensive alternative to the null-field 
condition, we propose to work with directive surface sources, which are obtained by an 
impedance boundary condition linking electric and magnetic surface currents densities and 
by shifting the equivalent surface current densities into complex space. These directive 
sources are integrated into the fast irregular antenna field transformation algorithm which 
achieves low numerical solution complexity by utilizing multilevel fast multipole like 
principles. Field transformation results based on real measurement data are investigated to 
illustrate the performance of the directive surface sources. 


Keywords: antenna measurement, near-field measurement, field transformation, antenna 
diagnostic, integral equation, fast multipole method 
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Outline TUTI 


> Introduction: Near-Field Antenna Measurements and Transformations 

> The field transformation problem with equivalent surface sources 

> Variations of the Huygens principle 

> Directive source modeling: Huygens radiators and sources in complex space 
> Results with measured and synthetically generated near-field data 

> Gaussian beam based plane wave translation operators 


> Conclusions 
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Dielectric Probe 
Waveguide 


Source: C. Koenen, et. al.: A volumetric Near-Field Scanner for Millimeter-Wave Antenna Measurements 
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Portable Antenna Measurement Systems (PAMS) 


a project with the European Space Agency (ESA) 
field transformation algorithms by TUM HFT 
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Automobile Antenna Measurements above PEC Ground TUT] 


> Cooperation with Antenna Technology Center (Europe) ATC GmbH 
= Hemispherical near-field antenna measurement facility 
= Metallic ground plane 


lS 





ANTENNA TECHNOLOGY CENTER 


http://www. antenna-tc.net 





Chair of High-Frequency Engineering Eibert November 2016 T 


Antenna Near-Field Measurements and Transformations TUT] 


YI = Measure AUT near-field for sufficient 
number of sample locations 


AUT 





= Representation of AUT by eguivalent 
>— by sources: 
o. e eigenmode expansions 
e equivalent currents 
e plane waves (far-field) 


representation by 


equivalent sources = Near-field samples related to 





probe equivalent sources (radiation operator) 
yo ON = Equivalent sources determined 
a f ay S b, by inverse problem solution 
/ \ 
= “2 by " Far-field/near-field obtained from 


\ 
I 
i equivalent sources 





Available transformation algorithms differ in 
their equivalent sources and radiation 
operator representations! 
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Why Inverse Problem Solution? — Influence of the Measurement Probe TUTI 
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= Real probe finite sized 
= Spatial/spectral filtering of measured field 
= Probe correction mandatory 

Filtering in spatial domain: 


UT) = |I| Poore F Fy) EF) 


VProbe 


Filtering in spectral domain: 
U (Fy) = Ep Werove(—k)- E(k ) dk? 


Wprobe(’) spatial probe weighting function 


W Probe (k ) spectral probe weighting function 


E(7) _ spatial near-field of AUT at 


oe a probe 
E(k) spectral near-field location 
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Radiation Model with Equivalent Surface Current Sources 


TUTI 


Consider the measured or simulated near-fields of an AUT or a scatterer (DUT). 


+ U (Tu) 
Virtual A O Tr 


Huygens f Serasi "" 


surface 17 } ee 
M 


y U (Tr) 


A 
wy. 7 
Huygens ii, 
surface `~ 
on object - = 
surface N , 


UED = [SÈ Worse F- 7u) |), 6; (EF) J, 0) + Gi F, F) M,C’) |da'dv 


VProbe 


Key question of this presentation: What are the most appropriate surface sources? 
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The Standard Form of the Huygens Principle TUTI 


> Electric and magnetic surface currents related to 
tangential fields: 


— — — 


i vA J,=nxH M, =-ñxE 


Virtual / a 
Huygens J " i > Mostly attributed to Schelkunoff/Love 
U de aes E f ~ 
a “se, M > Exclusion volume free of any field (null-field) 
X o e ; A A o 
Se A0 o 
da , 


Will we obtain this solution if we just solve the inverse problem as it is? 


S. Schelkunoff, “Some equivalence theorems of electromagnetics and their application to radiation problems," Bell System Technical Journal, vol. 15, no. 1, pp. 92—112, 1936. 
A. E. H. Love, “The integration of the equations of propagation of electric waves,” Philosophical Transactions of the Royal Society A, vol. 197, no. 287—289, pp. 1-43, 1901. 
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Equivalent Surface Current Discretization TUTI 


> Electric and magnetic surface current densities on triangular meshes 


— 


P, : vector basis functions 





RWG: Low order 
1st order Quasi Curl-Free 


> Derive li ion m by Method of 
Derive linear equation system by Method o 3 basis functions per triangle 


Moments procedure and solve it. 





p, =nx (AVA, E AVA, ) 
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Near-Field Measurement — HF907 by Rohde & Schwarz ILITI 
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spherical near-field measurement at 18 GHz measurement distance 2.73 m 


performed in the anechoic chamber at TUM 
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Transformed Far-Fields — HF907 ILITI 


Crosspol --=- FIAFTA curr. 
E-plane —— FIAFTA sph. 
—— NSI 


Diff. sph. 


---— FIAFTA curr. 
—— FIAFTA sph. 


— NSI 
Diff. sph. 


Relative far-field pattern in dB 
Relative far-field pattern in dB 





150 -100 -50 0 50 100 150 -150 -100 -50 0 50 100 150 
Y in degree Q in degree 


FIAF TA: Fast irregular antenna field transformation algorithm 
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Surface Current Densities J/M Without Further Condition — HF907 TUT] 


z undesired 


currents 








Chair of High-Frequency Engineering Eibert November 2016 15 


Electric Near-Field Due to -Hr9o7 ILITI 


Electric field magnitude in dBV/m Electric field magnitude ın dB V/m 
40 0.3 40 
30 
20 
z 10 
= 
À 0 





zın m 





xz-plane yz-plane 


Chair of High-Frequency Engineering Eibert November 2016 16 


Surface Current Densities J/M with Love Condition — HF907 TUT] 


J, in Am f Real J] in A/m = Imag J, in A/m_ 
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Variations of the Huygens Principle TUTI 


What kinds of surface current densities allow > Standard form, Schelkunoff/Love 


= 


correct representation of the radiation fields? J =nxH M,=-NxE 


e explicit enforcement required 
(Surface integral equation) 
e large computational effort 


i.e are suitable for the Huygens principle... 


+ U (ru) 


> J,,M, without further condition 


n vA 
Nase" = YI e works fine but redundant 

Virtual ,” eo : 
Huygens," 3 ` > Eitherof J, or M, alone 
surface _ 3 

+ Te d M, . e works but “bad operator” 

~ e nm ; : E E 
J ; >» Linear combination of J, and M, 
A 


e.g., M,=Z,(ixJ,) : imped. bound. cond. 


e approximate Love condition 
e so-called Huygens radiator 


S. Schelkunoff, “Some equivalence theorems of electromagnetics and their application to radiation problems,” Bell System Technical Journal, vol. 15, no. 1, pp. 92—112, 1936. 


A. E. H. Love, “The integration of the equations of propagation of electric waves,” Philosophical Transactions of the Royal Society A, vol. 197, no. 287—289, pp. 1-43, 1901. 
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Variations of the Huygens Principle TUTI 
The goal is to have surface sources, which radiate predominantly to the outside! 
> Love condition perfect, but with large computational effort (surface integral equation) 


> Impedance boundary condition M 4 =Z (x J a) gives directive basis functions 


— RWG Huygens radiator J, = LĒ, M, = I ,Z (ix B,) 


> Shifting the sources into complex space gives further directivity to the outside! 


e Huygens principle still valid 
e applicable to all types of sources as discussed before 
¢ well known within the framework of complex source beams/Gaussian beams 


i rss | ._—_ 
r =f a t JNA 


rea 


T. B. Hansen and G. Kaiser, “Generalized Huygens principle with pulsed-beam wavelets,” Journal of Physics A: Mathematical and Theoretical, vol. 42, no. 475403, pp. 1—33, 2009. 
T. B. Hansen and G. Kaiser, “Huygens’ principle for complex spheres,” IEEE Transactions on Antennas and Propagation, vol. 59, no. 10, pp. 3835—3847, 2011. 
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Radiation Fields of Rao-Wilton-Glisson-Radiators TUT] 
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Transformed Far-Fields — HF907 ILITI 


Source model: Electric and Huygens radiators with imaginary shift A = 0.242 
0 
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Electric Near-Field Due to Directive Sources — HF907 TUT 


Electric field magnitude in dBV/m Electric field magnitude in dBV/m 
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Surface Current Densities — HF907 TUT] 


J/M: no constraint, no imaginary shift Huygens radiators: A=0.2A 


> 1.50e-02 
< 1.50e-02 
< 1.39e-02 
< 1.27e-02 
< 1.16e-02 
< 1.05e-02 
< 9.33e-03 
< 8.20e-03 
< 7.07e-03 
< 5.94e-03 
< 4.80e-03 
< 3.67e-03 
< 2.54e-03 
< 1.40e-03 


< 2.70e-04 





> 8.80e-03 
< 8.80e-03 
< 8.14e-03 
< 7.476-03 
< 6.81e-03 
<6.15e-03 
< 5.48e-03 
< 4826-03 
< 4.16e-03 
< 3.49e-03 
< 2.83e-03 
<2.16e-03 
< 1.50e-03 
< 8.38e-04 
< 1.74e-04 


Max = 7.11e-02 Max = 3.22e-02 
Min = 2.70e-04 Min = 1.74e6-04 
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Electric Near-Field Due to Directive Sources — HF907 TUT 


Near-field along the AUT main beam axis dependent on the imaginary shift A 





electric currents Huygens radiators 
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The Spectral Domain Near-Field Operator TUTI 


UE) = SSS oo- E ED SfE EF IE) Gi EFM) jda 


VProbe 





—ų— a 


T,(k, X) : FMM translation operator 


U(r) — =j SP Prono (-k)- T, (k,F,) Sa) + Ma(k) (ak? : spectral representation 
It 
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The Spectral or Plane Wave Representation TUTI 


UT )=-j EGP Deok) = TE) Sak) +M a(k) ae? 
Ar 


LILI LIE Jalk)=(T—kk)- IE el) da! 


V Probe 


O : K O a 2 ik (7-7 ’ 
Probe radiation pattern M a(k) = iz— [Il x M , (7 eS da 
A 





AUT radiation pattern 
ik L | A - 
ZED UDP (kr PeP 


T, (k.F)= 
(KF) lee 





. Incident plane waves at 
field probe 


Plane waves 
radiated from AUT 
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Fast Multipole Near-Field and Far-Field Translation Operators TUTI 





FMM-Standard Translation Operator FMM-Far-Field Translation Operator 
= jkL ; : . pO! jkr|F 
T,(k,F) oe Pd (21+ Ih (kr) P,(k -?) i’ = F 0(9-3,)6(~-@) 
l=0 
> somehow directive, but with strong side lobes > very directive 
> all plane waves need to be considered > just one plane wave to be considered 


> applicable in far-field distance 


> good accuracy hard to achieve 


Field transformation problems have relatively large translation distances due to separation of 
sources and observation locations! 


This is in contrast to MLFMM solution of radiation and scattering integral equation solutions. 
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Fast Multipole Near-Field and Far-Field Translation Operators 


Translation operator with less side lobes than standard operator? 


T. J. Cul, W. C. Chew, G. Chen, and J. Song, 
“Efficient MLFMA, RPFMA, and FAFFA algorithms for EM scattering by very large structures,” 
IEEE Transactions on Antennas and Propagation, vol. 52, no. 3, pp. 759-770, 2004. 


R. J. Burkholder and D.-H. Kwon, 
“High-frequency asymptotic acceleration of the fast multipole method,” 
Radio Science, vol. 31, no. 5, pp. 1199-1206, 1996. 


e Ray propagation fast multipole method / windowed FMM 


Summation over the multipoles 7 with a windowing function 


e Fast far-field approximation fast multipole method 





eg" 
FF _ 
L = F 0(9-3,)6(9-Q) 
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Gaussian Beam Based Translation Operators 


Observation: Sources in complex space are directive. 
T. B. Hansen, 


“Translation operator based on Gaussian beams for the fast multipole method in three dimensions,” 
Wave Motion, vol. 50, pp. 940-954, 2013. 


T. B. Hansen, O. Bories 
“Gaussian translation operator in a multi-level scheme,” 
Radio Science, vol. 50, pp. 754-763, 2015. 


- jk| +d! 
ec 





za jk < a a ga 
x fpe 4 4 X\\ (7) (21 41) h® (kX) PB (kX) ak 
X +d tp tn |! |! 





T,(k,X) 
By choosing X+d =(X - jAÑ )+(d + jAX | 


a highly directive translation operator with very low side lobes can be constructed: 








- jk|X-+d| OL 
e jga Jk_kA(EX) a o | A A yn 
—— =dbe "+ e (-j)' (21 + Ih (kX — jkA)P(k-X)d'k 
X +d | tp 4r 2, 
T (k, X, 4A) 
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Gaussian Beam Based Translation Operators 
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0.05 


> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Gaussian Beam Based Translation Operators 
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> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Gaussian Beam Based Translation Operators 
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> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Gaussian Beam Based Translation Operators 
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> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Gaussian Beam Based Translation Operators 
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0.05 


> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Gaussian Beam Based Translation Operators 
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> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Gaussian Beam Based Translation Operators 
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> Considerably fewer plane wave samples needed 
> Extended translation ranges of benefit 

> Fewer MLFMM levels needed 

> Faster field transformation with less memory 

> Fully error controllable 


> Careful choice of parameters needed 
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Rohde & Schwarz HF907 Antenna — Transformation Performance TUT] 


> Mesh with 188 090 triangle 
> 393 452 Rao-Wilton-Glisson unknowns (one set) 


> 45 451 near-field sample points with two 
polarizations measured at every point 


> CPU: one core of Intel 17, 3.6 GHz 
> One iteration: 90s (standard FMM translator) 
12 s (Gaussian translator) 


7.5 times speed-up! 


> Typical iteration counts: 25 to 200 
(Dependent on the model and accuracy) 
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25 m (5002) Reflector Antenna — Method of Moments Solution TUTI 


59 329 410 RWG unknowns 
145 GByte RAM A 
6 GHz < 3.006-04 


<2./7e-04 





< 2? 54e-04 
<2.31e-04 
< 2.08e-04 
< 1.85e-04 
< 1.62e-04 
< 1.38e-04 
< 1.156-04 
2 9 23a-05 
< 6.92e-05 
< 4.622-05 
<2.31e-05 
< 0.00e+00 


Max = 7.81e-04 
Min = 8.832-08 





E 4-Element Yagi Feed f 
d with delta gap voltage source 


L NS A, 
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25 m (500A) Reflector Antenna — Synthetic Near-Field Data TUTI 


regular sampling, z = 7 m (1 m above aperture) 

4/2 sample distance | | | 

5 764 801 measurement locations Normalized IEI in dB in scan plane 
(two orthogonal polarizations each) obtained from MoM solution 


0.25 





500 1000 1500 2000 
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25 m (500A) Reflector Antenna — Transformation Results TUTI 
Regularly sampled near-field data 


Equivalent surface currents from FIAF TA transformation 
z = 6.13 m (slightly above aperture) Copol, E-plane 


MoM ref 
FIAFTA 





Electric far-field in dB V/m 


Error in dB 


Max = 1.7 78-03 
Min = 7.46e-08 
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e 41 775 501 surface current unknowns 


e CPU time on one core of Intel i7, 3.6 GHz: about 1378 sec per iteration (99 iterations, final result) 
e RAM about 106 GByte (electric currents) 
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25 m (5002) Reflector Antenna — Synthetic Near-Field Data TUTI 


strongly irregular sampling, z= 7 m (1 m above aperture) 
A/2 sample distance on average 
5 /64 801 measurement locations 

(two orthogonal polarizations each) 


-19.3 


yın m 





-19.5 


-30 -29 -28 -27 -26 -25 -24 
X IN M 





-0.9 -0.7 xin m -0.3 -0.1 
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25 m (500A) Reflector Antenna — Transformation Results TUT] 


Copol, E-plane 
—— Regular 


Relative far-field pattern in dB 





“80 60. 40 20 “0 20 40 60 80 
Q in degree 





Y 


4-Element Yagi Feed 


with delta gap voltage source 
Directive sources very helpful for 
irregular sampling on order to achieve 
x "O low error! 
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Conclusions TUT] 


Near-field transformations with directive sources: 


> 


> 


Two ways: Huygens radiators and imaginary shift into complex space 
Huygens principle still fully applicable 

Improved diagnostic capabilities 

Improved conditioning of resulting equation systems 


Improved transformation performance with Gaussian beam translators 


Chair of High-Frequency Engineering Eibert November 2016 43 


References TUTI 


A. D. Yaghjian, An Overview of Near-Field Antenna Measurements, IEEE Transactions on Antennas and Propagation, vol. 34, no. 1, pp. 30—45, 1986. 


C. H. Schmidt, M. M. Leibfritz, T. F. Eibert, Fully Probe-Corrected Near-Field Far-Field Transformation Employing Plane Wave Expansion and Diagonal 
Translation Operators, IEEE Transactions on Antennas and Propagation, Vol. 56, No. 3, pp. 737—746, March 2008. 


T. F. Eibert, C. H. Schmidt, Multilevel Fast Multipole Accelerated Inverse Equivalent Current Method Employing Rao-Wilton-Glisson Discretization of Electric 
and Magnetic Surface Currents, IEEE Transactions on Antennas and Propagation, Vol. 57, No. 4, pp. 1178—1185, April 2009. 


C. H. Schmidt, T. F. Eibert, Multilevel Plane Wave Based Near-Field Far-Field Transformation for Electrically Large Antennas in Free-Space or Above 
Material Halfspace, IEEE Transactions on Antennas and Propagation, Vol. 57, No. 5, pp. 1382—1390, May 2009. 


T. F. Eibert, Ismatullah, E. Kaliyaperumal, C. H. Schmidt, Inverse Equivalent Surface Current Method with Hierarchical Higher Order Basis Functions, Full 
Probe Correction and Multilevel Fast Multipole Acceleration (Invited Paper), Progress in Electromagnetics Research (PIER), Vol. 106, pp. 377—394, 2010. 


T. F. Eibert, Emre Kilic, C. Lopez, R. A.M. Mauermayer, O. Neitz, G. Schnattinger, Electromagnetic Field Transformations for Measurements and 
Simulations (Invited Paper), Progress in Electromagnetics Research (PIER), Vol. 151, pp. 127—150, 2015. 


S. Schelkunoff, “Some equivalence theorems of electromagnetics and their application to radiation problems,” Bell System Technical Journal, vol. 15, no. 1, 
pp. 92—112, 1936. 


A. E. H. Love, “The integration of the equations of propagation of electric waves,” Philosophical Transactions of the Royal Society A, vol. 197, no. 287—289, 
pp. 1—43, 1901. 


T. B. Hansen, G. Kaiser, Generalized Huygens Principle with Pulsed-Beam Wavelets,” Journal of Physics A: Mathematical and Theoretical, vol. 42, no. 
475403, pp. 1—33, 2009. 


T. B. Hansen and G. Kaiser, Huygens’ principle for complex spheres, IEEE Transactions on Antennas and Propagation, vol. 59, no. 10, pp. 3835—3847, 
2011. 


T. F. Eibert, D. Vojvodic, T. B. Hansen, Fast Inverse Equivalent Source Solutions With Directive Sources, IEEE Transactions on Antennas and Propagation, 
Vol. 64, No. 11, November 2016. 


T. B. Hansen, Translation Operator Based on Gaussian Beams for the Fast Multipole Method in Three Dimensions, Wave Motion, vol. 50, pp. 940-954, 
2013. 


T. B. Hansen, O. Bories, Gaussian Translation Operator in a Multi-Level Scheme, Radio Science, vol. 50, pp. 754-763, 2015. 


Chair of High-Frequency Engineering Eibert November 2016 44 


